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Abstract 

Michor and Mumford have shown that the distances between planar 
curves in the simplest metric (not involving derivatives) are identically 
zero. We consider two conformally equivalent metrics for which the dis- 
tances between curves are nontrivial. We show that in the case of the 
simpler of the two metrics, the only minimal geodesies are those corre- 
sponding to curve evolution in which the points of the curve move with 
the same normal speed. An equation for the geodesies and a formula for 
the sectional curvature are derived; a necessary and sufficient condition 
for the sectional curvature to be bounded is given. 

1 Introduction 

The purpose of this paper is to study the most basic properties of some of the 
simplest Riemannian metrics suggested by applications to Computer Vision. 
The problem is to understand and quantify similarities and differences between 
object shapes and their individual variations. At a fundamental level, the 
problem is to construct appropriate metrics on a space of closed surfaces in 
R 3 . A simpler version of the problem is the construction of Riemannian metrics 
on a space of closed planar curves. The choice of a metric depends on the type 
of similarity that is being considered. In their seminal paper [2], Michor and 
Mumford analyze two Riemannian metrics on a space of closed planar curves. 
Surprisingly, the Riemannian distance between any two curves in the simpler 
of the two metrics, an _ff°-metric, turns out to be zero. To remedy this, they 
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add a curvature term to the metric. Alternative construction of Ricmannian 
metrics on the space of curves of constant length and parametrized by cither 
the orientation of the tangent vector or the curvature is given by Klasssen 
et al in [3]. A method for deforming one curve onto another by minimizing 
an approximate Hausdorff distance between them is described in [1]. Below, 
we analyse two conformal variants of the 7?°-metric of Michor and Mumford. 
In §3, we derive upper and lower bounds for distances between curves and 
show that these metrics behave like L 1 metrics. A key issue is the existence of 
minimal geodesies. We show that for the simpler of the two conformal metrics, 
the minimal geodesies correspond exactly to those curve deformations in which 
the points of the curve move with the same normal speed. In the case of the 
second metric, no geodesic is minimal if the length of the curve is less than a 
certain threshold; the question of minimality when the length of the deforming 
curve is equal or greater than threshold is still open. We provide a partial 
answer in the form of a necessary and sufficient condition for the boundedness 
of the sectional curvature in §5. An equation for the geodesies is derived in §4. 

This paper began with an analysis of the simpler of the two conformal 
metrics considered here. As it was being written, the author became aware 
of the work of A. Yezzi and A. Mennucci [4,5] in which they have proposed a 
more general formulation. Since our analysis applies to their formulation as 
well, we have included their formulation in the analysis below. 

2 The Framework 

The basic space considered by Michor and Mumford is the orbit space 



of the space of all C°° embeddings of S in the plane, under the action by 
composition from the right by diffeomorphisms of the unit circle. It is contained 
in the bigger space of immersions modulo diffeomorphisms: 



Let 7r : /mmfS 11 ,! 2 ) — ► i?i(S' 1 ,R 2 ) be the canonical projection. The simpler 
of the two metrics considered in [2] is an iJ°-metric defined on Imm(S 1 ,R 2 ): 



where c : S 1 — > R 2 is an immersion, defining a point in Imm(S 1 , R 2 ), h,k 
£ C°°(S 1 , R 2 ) are the vector fields along the image curve, defining two tangent 
vectors in Imm(S 1 ,R 2 ) at c, and eg — dc/d9. (h-k) is the usual dot product in 
R 2 . Sometimes for the sake of clarity, we will use the notation a ■ b even when 
a, b are scalars. Let n c denote the unit normal field along c. If we identify R 2 
with the complex plane C, then, n c — ice/ \ c e\- For any C Q , C\ G Bi, consider 



BeOS' 1 , R 2 ) = Emb(S 1 7 R 2 )/Diff(S 1 ) 



Bi(S\R 2 ) = Imm{S\R 2 )/Diff{S 1 ) 




(1) 
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all liftings c ,ci to Imm(S 1 ,M 2 ) and all smooth paths t i— ► (9 i— » c(t,9)), 
< t < 1, in Jmm(S 1 ,l 2 ) with c(0, •) = c Q and c(l, •) = c\. Let c t denote 
dc/dt and q 1 = (c t • n c ) n c . The arc-length of such a path c is given by 

/ y/G% (c t ,c t )dt 
Jo 

Michor and Mumford show that for any two curves in B i (S 1 ,M. 2 ) , 

dist G o (Ci , C 2 ) = def inf J yj G° (ct,c±-)dt = 

and strengthen G° by defining 

G*(h,k)= [ (l + An 2 c )(h-k)\c e \d9 
Js 1 

where n c is the curvature, defined by the equation (jf^y) = * K cCe = k c \c$\ n c . 

An alternative is to consider conformal transformations of G° . They have 
the form 

G*(h,k) = *(c) / (h-k) \c \d9 
Js 1 

where $ is a £>i//(S' 1 )-invariant function on Imm(S 1 ,R 2 ). For example, 

$(c) = / \ce\d0 or $(c) = / (l + Ak 2 ) | Ce |cf0 

Based on considerations of stability, Yezzi and Mennucci have proposed a con- 
formal factor of the form e Al where i = / sl \cg\d9, or more generally, a function 
ip(£) of I. In this paper, we consider only these conformal factors: 

G%(h,k)=<p(e) f {h-k)\c e \d9 (2) 
Js 1 

If c(t, 9) is a smooth path in Imm(S 1 , M 2 ) connecting c Q , ci, let 

L G «{c) = j^Gt (ct-ct)dt (3) 

If C ,Ci e B^S 1 ,®?), define disi G ^ (Ci,C 2 ) as follows. Consider their lifts 
c D ,ci to Imm(S 1 , R 2 ) and all paths c(i, 0) such that 7r(c ) = C and 7r(ci) = 
Ci. Then 

disi G *> (Ci, C 2 ) - inf L G „(c) (4) 
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3 Bounds on distcv 

If c is a path connecting curves C, E, let a(c) denote the area swept out by c 
in R 2 . For a path c(t, •), let £ max (c) = max. t £(c(t,-)). The following theorem 
characterizes the L 1 -type behavior of the metrics G v . 

Theorem 1 // ip{£) = i, then, 

dist Gv (C,E) = in£a(c) (5) 

C 

If <p(£) = e M , then, 

inf y/Aea(c) < dist G * (C, E) < inf y/Tee M ^ {c)/2 a{c) (6) 

c c 

We first prove a series of lemmas. 
Lemma 2 

dlStGV {C ' E) ~ { inf c VTeaic) if ^ = e Ai 
Proof: For any path c, 



Js 1 



dt 



> 



> 



o 



\supp{cf)\ 



\ct\\<»\ 



dt 



mm — 



> 



f a(c) 
1 V^fc 



Six [0,1] 



det dc(t,6)\d6dt 



(c)_if = £ 
Aia(c) if = e A£ 



Q.E.D. 

David Mumford observed from the formula for the sectional curvature that 
the geodesies along which \supp(c^~)\ < £ if ip(£) = £ and < 1/A if ip(£) = e M 
may not be minimal. Such a possibility can be heuristically seen from the 
inequality 



L g >p{c) > 



V (£) 



\supp(cf)\ 



c e \dO 



dt 



which suggests that while traversing a given area, one should try to minimize 
ip(£) /\supp(ct )|. The key point is that we can increase \supp(c^)\ indefinitely 
by replacing the part of the curve supporting q 1 by a saw-tooth shaped curve 
of high frequency and small amplitude. When ip(£) = £ and \supp{ q*-)| < £, 
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we can increase \supp(c^)\ so that £/\supp(c^~)\ tends to 1. When ip(£) = e M 
and \supp(c^)\ < 1/A, we can force e A<[SupplyC t^ /\supp(c^)\ to equal its unique 
minimum Ae by making \supp(c^)\ equal 1/A. (In the case of the metric G° , 
<p(£) = 1 so that <p(£) /\supp(c^)\ tends to 0.) In order to obtain an upper 
bound for a general path, we break it up into a series of tiny bumps. When 
ip{£) = £, this method gives an upper bound for distct (C, E) which coincides 
with the lower bound. When tp(£) = e M , the larger the value of \supp(c^~)\, 
the greater the divergence between the upper bound obtained by this method 
and the lower bound since it is more efficient to create a large bump all at once 
instead a series of tiny bumps. 
Rectangular Bumps 

Let c a : S 1 — ► R 2 be a smooth and free immersion. Let C a be the cor- 
responding curve in R 2 . Let c Q be parametrized by the arclength so that 9 
parametrizes the scaled circle S\ where £ Q is the length of C Q . For any func- 
tion u(9), let u' denote du/d9. Let n D denote the normal vector ic' Q . Let k 
denote the curvature of c . 

Fix small positive numbers 5 and e such that 5 < £ and e H^Hoo [05] << 1- 
Construct a "rectangular" bump, over C as follows: 

C c (9) + en H0<9<5 
ci(0)=< {c (9) + sn o \0 < s < e} if 9 = 0,6 
[ c (9) otherwise 

Let C\ be the corresponding curve in M 2 . The following lemma is inspired 
by a comment of David Mumford that the Michor-Mumford "teeth" construc- 
tion [2] may be used to show that the obvious path for creating such a bump 
is not minimal. 



Lemma 3 (i) If <p(£) = £, 



distG*>(C ,Ci) < 



1 + e|M|oc.[o,5- 1 " 



.1 e ll K o||oo : [o,(5] 

(ii) If cp(£) = e M and 6 < 1/A, 



(area of the bump) 



dist G v(C ,Ci) 



< 



1 + e||«o||oo,[o,5] 



1 - e|k | 



[o,5] 



3/2 



l+«ll«oll 



e A(e o +2e-6)/2y Ae i-^oU^ [o , a] ( bump Qrea 



Proof: We prove the lemma using a modification of the "teeth" construc- 
tion of Michor and Mumford [2]. If ip(£) — £, choose A < 1/5. Approximate 
C\ by a "trapezoidal" bump C\ as follows. Replace Co in the interval [0,(5] 
by a saw-tooth curve of height -q and period ^ such that its length equals ^ . 
This is done by growing teeth on Co in time r\. Move the saw-tooth curve at 
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unit speed along the normals no keeping its end-points fixed, until it touches 
the upper edge of the bump. Finally, retract the teeth in time rj. Formally, 
define a path c(t, 9) — c o (0) + f(t, 9)n Q where f(t, 9) is defined as follows. 



f(t, 6) = 0, < t < e and 5 < 
For < t < rj, < k < m - 1, 



/(*,*) = 

For r] < t < e — r], 

' e 

f(t,6) = < 
For e — rj <t < e, 
f(t,6) 



t {2k + 2-^f. 



t{t-rt)+r)(e-r)-t) 2m0 
e-2r/ ' S 

e(t—T))+ri(e-ri-t) 



2k < 8 < 
2m — S — 2m 



2fc+l 



2fc+l 
2m 



< ? < 



2fc+2 
2m 



o < 

1 

2m 



< 



1 

2m 

< 1 



1 

2m 



e-2rj 



2m (1 



2m 



< ? < 1 



2me9 

elt-(e- n )] + (e~t)f(e-r h e) 

2me (1 - f 



0< | < ^ 

j_ < 4 < i _ j_ 

2m — o — 2m 



2m 



< 



< 1 



c =c' + f'n - fK c' = (1 - fn )c + f'n 

\c'\ = VTi - M + f' 2 
_ -f'cfp + (1 - fKo)n 
Ic'l 



Let 



<h = f t n 
cf- = (c t ■ n) 



2 (1 - ]n Q ff t 







1 < < 



' (l + e||« || 00 , M] ) 2 + /' 2 
(l-e||Ko|]oo,M]) 2 + /' 2 

1 + e| \k I |oo, [o, 5] 



1 - e||«o||oo,[o,5] 



Choose m and rj such that \c'(rj, 9)\d9 = -j. Note that as m — > oo, 
r? -> 0. 
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Estimates when < t < r\: Since \ is independent of 9 and |/| < e, 

/(HWU^F + f < \c'(v,o)\ < ^(i + e\K\U[o, s] ) 2 + f /2 

5^/(1 -e||K ||oo,M]) 2 + / /2 < f l c '(^ 0)\d9 = 1 < ^(l + 6|KIU M] ) 2 + .f 2 
We also have 

/(HWUiW < 0)| < ^(l + e||«o||oo, M ]) 2 + /' 2 

+ Jmumi) 2 + / ,2 < < - eii«oiioo, M] ) 2 + f 2 



H^ |c ' M)l ^ 



1(c) = j ° \c'\d6 

e AUc) < e A{toS) e D 



Since |/ t | < 1, 



j o °\ci\ 2 \c'\de = ^ (1 (j | o)2/t2 rfe<(i + e ||« ||oo,M]) 2 /3^ 



lim 

m— »oo 



-,1/2 



^(£( C )) / l^lVMtf 



«it = 



Estimates when r\ < t < e — r\: Estimate for \c'(t, 9)\ is the same as in the 
interval [^,8(1 — 5m)] since the curve has the same shape. In the intervals 
[0,£] and [6(1 -^),S], ^ < \f(t,6)\ < Therefore, \c'( V ,8)\ < 

\c'(t,6)\ <1 + e| Klloo,[o,a] + ^- 



*(c) = ^ ° \c'\d6 < (£ a -6)+(l + e||« ||oo,[ ,5] + 

1 + e||Ko||oo,[o,5] 1 



2me\ S f3 
5 ) m A 



lim (.(c) < (£ - 5) + 2e + 



1 - e|k | 



[o,S] 



1 +«ll^°lloo,[o,a] 

lim e Ae(c) < e ^+ & - s ) e I=;,ra tM 

m— »oo 

We also have < ^f^- Therefore, 



jf /0 |^| 2 |c'|^ < (1 + e||« |U, M] ) 2 ^ 
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lim f " <p(l{c)) f°\ci\y\d6 

wo % Jo 



1/2 



dt 



1 + e||/c 



l + e||Ko||oo,[ ,j] 



< < 



if ip(£) = i 



v ° ; l-e||Ko||oo,[o,«] 



1/2 



I ! I £||Ko||ocJo,d]j y l-e||Ko|| 00 , [o „ 5] e 



ll"oll 



Estimates when e — r] < t < e: The path c(i, 0) in the interval [e — r), e] is 
essentially the same as that in [0, 77] and 



lim 

m — >oo 



-,1/2 



|^lVl<» 



dt = 



Since 



(area of the bump) = I I (1 — tn o )d0dt > (l — e||«o||oo,[ ,5]) ^ e 

the lemma is proved in the case when yj(^) = e A£ . It is proved in the case when 
ip(£) =£by letting A -> 0. Q.E.D. 

In Lemma 3, we may replace the single bump by a finite number of disjoint 
bumps of height e and total length 5. The proof remains unchanged except 
that we must replace 2e in the formula by 2ke if k is the number of bumps. 
The function / in each individual bump may be positive or negative. 

We prove the theorem by approximating the path by a series of small 
rectangular bumps. The error of approximation may be made arbitrarily small 
by the following lemma. Define a Frechet metric on Imm(S 1 , R 2 ) and 
B^S* 1 , R 2 ) as follows. If c ,ci are points in Imm(S 1 ,R 2 ), let 

doo (c , ci) = sup I c (0) - ci (6) I 

If Co, Ci e R 2 ), let 

d 00 (C 1 ,C 2 )= inf rfcx>(c ,ci) 

{c ,Ci|7r(c )=C' ,7r(ci)=C'i} 

Lemma 4 For any pair 

C U C 2 e Bi^.M 2 ), dist G * (Ci.Ca) < ^(Ci, C 2 ) • max{^i), ^ 2 )} (7) 
w/iere ^ = £(Ci), i = 1,2. 
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Proof: Letci c 2 be lifts of Ci, C 2 to 7rom(S\R 2 ). Let c(t, 9) = (1 - t) Ci (0)+ 
tc2 (9) be a path connecting them. Then, \cg(t)\ < (1 — t) \c\ t e\ + t\ c 2,e\ and 
hence, £{c{t)) < max{^(Ci), £(C 2 )}. Moreover, c t = c 2 — Ci. Therefore, 

dista-c (Ci,C 2 ) < infi Gv (c) 

c 

< inf (sup|ci(0) - 02(^)1 }max{^i),^ 2 )} 

{pairs ci,c 2 } J 

Q.E.D. 

The polygonal approximations used in the proof of the theorem lie on the 
boundary of /mmfS 1 ,! 2 ) and B^S 1 , M 2 ), and Lemma 4 extends to them. 
Proof of the theorem: 

Consider a path c (i, 6>) connecting C and Since the absolute curvature 
of the curves c(t, •) is uniformly bounded by a constant K, each curve has a 
tubular neighborhood of width which is bounded from below. Choose e and 
a sequence 

= t < h < ■ ■ ■ < t N -i < t N = 1 
such that eK « 1 and, for < k < N, c (tk+i, •) is in a local chart of c (tk, ■)■ 

c(tfc+i,0) = c(tk,9) + f k {6)n k 

where c(tfc,-) is parametrized by the arclcngth, nk is the normal vector field 
of c(t fe ,-) and l/fel < e. Let F = max{||/^ (0)^ |0 < fe < N}. Let C fc denote 

7r(c(t fc ,-)). Let l k =t(C k ). Letl k = (l + eK + F)£ k . 
Choose S such that 

max \Jt k ip{i k ) ■ FS < 

We now estimate the distances distc v {C k , Cfc+i). Consider the path seg- 
ment [c D , ci] from c to c\ in the local chart at c . Divide the range of 9 
into intervals of length 8 . Replace f a by a piecewise constant function f a 
whose value in each subintcrval equals the average of f a over that interval. Let 
C be the curve defined by f . The Frechet distance between C\ and C is 
< FS. The sum of the jumps in /„ is < Fl Q . Since \c' D \ = |1 — ck \ < 1 + elf, 
^Co) < (1 + eK + F)£ = £ . Therefore, 

dist Gv (d,C ) < max{i <p(io),iMii)} ■ FS < 

Let a([c ,ci\) denote the area swept out by the path c during [0, £i] . The 
area between C and C equals the area between C and C\ which in turn is 
less than or equal to a([c ,ci]). The curve C consists of a series of bumps 
over C - Traverse the bumps sequentially, taking care to retract the common 
edge of each bump with the previous bump before going to the next bump. 
Retracting a common edge can be done without incurring any cost. 
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By Lemma 3, if <p(£) = £, 

dist Gv , (C , C ) < ( YS^k ) a ^ Co ' 



1 + eif N 2 



and if ip[t) = e Ae , 

3/2 



dist G , (C , C ) < e A ^y^Ae^a([c 0l Cl ]) 

dist G * (C ,Ci) < dist G ~p (C ,C) + 

Similar estimates hold for dist G v (Cfc, Cfe+i) for < k < N. 
Therefore, if ip{£) = £, 



dist Gv (C, E) < y 1 _ eK J a ( c ) + e 
and if tp(£) = e M , 

dist Gv (C,E) < ( e M ™ (c)/2 \pte£&a(c) + t 



1-eK 

Since e is arbitrary, we have 



dist Gv [C,E) < \ , ;/U , ! , ux 2 ,-r- rl , . ir , w 



a(c) if = £ 

/ Iea(c) if <p(£) = e 1 

Q.E.D. 

For any oriented curve C or , define the integer- valued measurable function 
w G on R 2 by: 

w G (x,y) — the winding number of C around (x,y) 

and let 

d\Cl r \C° 2 r ) = f \w Cl - w C2 \dxdy 
It is shown in [2] that for any two oriented curves C^ r , C% r , 
d b (Cr,CD< „ rt min. a(c) 

all paths c joining G1.G2 

Therefore, we have 

Ofe S i G „ (Gi ,G 2 ) > I Jfe# {C cr^ C o r) if ^ = e ^ I 8 ) 
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Corollary 5 (Existence of minimal geodesies) If <p(£) = £, then the only min- 
imal geodesies are the paths along which \c^\£ is constant. 



Proof: Since the inequality 



L Gv (c) > 



\ct\\c e \de 



dt = a{c) 



is an equality if and only if | does not depend on 9 , that is, ^ \cj~ | =0 (the 
case of " grassfirc" ) , if c(t) is a minimal geodesic, \c^\ must be independent of 
9. Let c(t,9) be a path connecting Ci,C2 such that \c^\ is independent of 9. 
After reparamctrization if necessary, we may assume that CfCe = . Following 
[2], we let w(c) be the 2-current defined by the path c(t, 9). Since c(i, 9) is an 
immersion, 

d\Cl r ,C° 2 r ) = [ \w(c)\dxdy= [ | det dc(t, 9)\dBdt = a(c) 

Jm 2 is 1 x [0,1] 

Therefore, Lg* (c) is the minimal distance between C\,Ci. For c(i, 9) to be a 
geodesic path, reparametrize t such that the infinitesimal arc-length \c^~\£ is 
constant along the path. Q.E.D. 

Corollary 6 Suppose <fi(£) = e Ai and c(t,9),0 < t < 1, is a path connecting 
Ci, Ci. Assume that \c^\ does not depend on 9 and £{t) < 1/A for all t. Then, 



dist G ACi r ,C2 r ) = ^Ted\Cl\Cl r ) 

It follows that the path c(t, 9) is not minimal. 

Proof: Break up the interval [0, 1] into small segments of length e and 
apply Lemma 3 with 5 = £ . (The proof of the lemma extends to this case 
after minor modifications.) Calculate the length of the new path c £ applying 
the construction of Lemma 3 to each of the segments.. We get lim £ ^ L GV (c £ ) = 
\/Aea(c) as in the proof of Theorem 2. On the other hand, 



Lg*(c) = J y \ci\\c e \d9^ dt = J ^f^±\ci\£dt 



> \/min ^p-a(c) = \ /f Aea(c) 

Q.E.D. 



4 Geodesic Equations 

We reproduce calculations in [2] mutatis mutandis. Let 1 1 — > c(t, •) be a path 
in Imm(S\R 2 ). 
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Equation of Geodesic in Imm(S 



1 TB>2\ 



W)\ce\ct) t 



\c t \ Cg 
M 



(9) 



where ip' denotes d<p/d£. 

Proof: We calculate the first variation of the energy of the path to obtain 
the geodesies:. 



E G * (c) 



1 



b r 



ip(e)[ l ( Ct -c t )\cg\de 

S 1 



dt 



d s \ ip{i) = ip'{£)d s \ / (Cg-Cg)*d9 

= / (c s ■ 

J si Vl c e I 

= -<p'(£) / (c s • in c cg) dO 
Js 1 



ds\o\ I {c t -c t )\ce\de= [ ( Cst -c t )\cg\d8 + l [ \c t \ 2 ^'^ dO 
2 Js 1 Js 1 2 Js 1 \ c s\ 



[ (c st ■ Ct) \cg\d6 ~ \ f {c s 
Js 1 1 Js 1 



\c t \ 2 Cg 
M 



)d8 



d s \ E Gv (c) = ~<p'(e)J a {c s -iK c c e )de^ Qj^ |c t | 2 \c e \d6^ dt 



u,i^M )ds)dt 



+ la ^ [j sl ^- C ^\ C0 \ d6 ) dt 

The last term = / (c s • (p(£)c t \c e \ t ) dOdt 

Ja Js 1 



12 



Therefore, 



d s \ E Gv (c) = - / / (c s -F)d9dt 

Ja JS 1 

where F = (<p{£) \c e \c t ) t + (J \ Ct \ 2 Mc ^ 

i (\ot\ 2 c e \ 
2 ^ \c e \ ) g 

Q.E.D. 

Equation of Geodesic in B i (S 1 ,R 2 ) : To obtain the geodesies in B i (S 1 ,M 2 ), 
we may assume that (q • eg) — 0. Write Ct as iacg/ \cg\. After substituting this 
in Eq. (9), split the equation into its component along eg and n c . The former 
vanishes identically. The normal component of the left-hand-side is 

(ip(£)\cg\a) t =(p(£)a t \cg\-(p(£)K, c a 2 \cg\~ip'(£)a\cg\ / an c \cg\d9 

Js 1 

where we have used the formulae \cg\ t — —an c \cg\ and £ t (c) = — J gl an c \cg\ d6. 
The normal component of the right hand side is 

-\<p{£)n c a 2 \c e \ - \y'{t)K c (^J a 2 \cg\ d9^j \c g \ 
Therefore, the equation of the geodesic in B i (S 1 ,R 2 ) is 



f - (a 2 -a 2 ) + a ■ ai^ if ip(£) = £ 

% (a 2 - (A£)a?^ + (A£)a • a7^ if <p(£) = e M 



(10) 



where for any function f(t, 9), / denotes the average j f sl f \cg\ d6. 

As an example, consider the case of concentric circles [2], c(t, 9) = r{t)e ie , r a - 
r(0),ri = r(l) We have k c — 1/r and a = —r t . Substituting these in Eq. 
(10) when ip(£) = £, we get —ru — r 2 /r or (r 2 ) u = 0. Therefore, r 2 (t) = 
tr\ + (1 — t)r 2 . This example is a special case of the curve evolution by "grass- 
fire" in which a is independent of 9. We have a 2 = a 2 and the equation of the 
geodesic reduces to a t — a-an^ — —a£ t /£ and hence (a£) t = 0. Therefore, a£ = 
a constant. By substituting in the equation for the length of the geodesic, we 
find that a£ = the length of the geodesic. When ip(£) = e Al 1 the equation of 
the geodesic in the case of concentric circles is (r 2 ) u = rf(l — 2irrA) which is 
zero when the perimeter of the circle equals 1/A, marking the unique inflection 
point of the function r 2 (t) . 
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5 Sectional Curvature 

A formula for the sectional curvature may be derived exactly as in [2] by 
means of local charts. Let c : S 1 — > M 2 be a smooth and free immersion, 
an element in the space of free immersions, /mm/fS 1 ,! 2 ). Bij(S , R ) = 
7T {lmmf(S 1 , R 2 )) . Let C = 7r(c). As before, let c be parametrized by the 
arclength so that 8 parametrizes the scaled circle S} where I is the length of 
c. Let 

V> : C°°(Sj, (-€,e)) -» 7mm / (S* 1 ,R 2 ) 
iKf){6)=c{O) + f{0)n c {6) 
tto^j: C°°(Sl (-e, e)) -> B ii/ (S 1 ,R 2 ) 

be a local chart centered at C. For any function u , let v! denote its derivative 
so that <p' — dip/dl and if u is a function on S\ , u' = du/d9; . We have the 
following formulae from [2]: 

= (1 - fn c )c' + f'n c 
(1 - f Kc )n c - fcJ 

nw) Vd-f- c ) 2 + f' 2 

For h € C°°(Sj,R), ft - n c € T w) Imm f {S l , R 2 ) and 
Let 

/■ (1 - fK c ) 2 hk ,„ 

<h,k> f = / - ; ' d9 

Js} y/(l - /K C ) 2 + /' 2 
= / /»fc(l-/ /6c -l/'2 + 0(/ 3 ))de 

*/ = / V(l - /^c) 2 + Z' 2 ^ = / (1 - /k c + ^/' 2 + O (f) )d6 
Js\ Js\ 1 

< h,k > will be denoted simply as < h, k > and \ \h\\ 2 —< h, h >. The metric 
at the point in Bij(S 1 , R 2 ) corresponding to / is given by 

G*(h,k) = <p(lj) <h 1 k> s 

We now outline the calculation of the Christoffcl Symbol and the sectional 
curvature. (For more details, see [2].) Differentiating the metric in direction 
j, we get 

d G){h,k) = tfitrfdjlf) <h,k> f +tp{l J )d j < h,k> f 
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- < j, k c > + < f',f > +0(f) = - < n c + f" ,j > +0(f) 



dj <h,k>f=-< hk,jK c + f'f > +0(f ) 

Therefore, 
d 3 G v f {h,k) 

= <p'{i f ) (- < j, K C > + < f',f >) <h,k>f- ifilf) < hk,j Kc + f'f > +0(f) 

= -<p'(£ f ) < K c + f",jxh,k> f -<p(£ f ) < hk,jK c + f'f > +0(/ 2 ) 
The Christoffel symbol satisfies the identity 

2G*(r f (h,k),j) = d j G ip f (h,k)- d h G*(k,j)- d k G v f (h,j) 
Therefore, 

2G£ (T f (h, k) ,j) = <f/(l f ) <k c + f", h > < k,j >f + <p(£ f ) < kj, hn c + f'ti > 
+ <p'{£f) < n c + f",k> < h,j > f +tp(l f ) < hj,kn c + f'k' > 
-if 1 {If) <n c + /", j ><h,k> f - <p(£ f ) < hk,jn c + f'f > +0(f 2 ) 

At the center / = 0, this simplifies to 

2G£ (r o (h, k) ,j)=< tp(£)K c hk+tp'(£) [< k c , h > k+ < k c , k > h - < h,k> k c ] , j > 
Since G% (T (h, k) , j) = tp(£) < T (h, k) , j >, we have 

r o (h, k) — -n c hk + f j^r [< k c , h > k+ < n c ,k > h— < h, k > k c ] 
2 2<fi{£) 

Next, we calculate the second derivative: 
d m djG^(h, k) 

= d m {<p'(£ f ) (- < j, Kc > + < f',f >) <h,k> f - <p(t f ) < hk,j Kc + f'f > +0(f 2 )) 
= <p"(£ f )(<m,K c > - < f',m' >)(<j,K c > - < f',f >) <h,k> f 
+ ip'(£f) < m',f ><h,k> f +ip'(£ f )(< j,K c > - < f',f >) < hk,mn c + f'm! > 
+ <p'(£ f )(< m,K c >-< f',m' >) < hkjKc + f'f > - <p{£ f ) < hk,m'f > +0 (/) 

At the center, 

d m d,Gt{h,k) 

= ip"(£) < to, k c >< j, k c >< h, k > 

+ <p'(£) (< m',f > < h,k > + < j,K c > < hk, mn c > + < m, k c > < hk,jn c >) 
- ip(£) < hk,m'f > 
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Let P(m, h) be a tangent plane at c, spanned by normal vector fields 
m(8), h (8) along C. Assume that m, h have been normalized so that 

<p(£) ||m|| 2 = ip(£) \\h\\ 2 = 1 and < to, h >= 

The sectional curvature at C is given by the formula 

k c (P(m, h)) = d m d h G%(h, to) - ^ {d m d m G£(h, h) + d h d h G%(m, to)] 

+ Gi (r (h, to) , r G (h, to)) - Gi (r G (to, to) , r (h, h)) 



d m d h G%{h, to) - - \d m d m G*{h, h) + d h d h G*(m, to)] 

= <p' (< m,K c > < mh, hn c > + < h, k c > < m/i, mn c >) — <^(f) < mh, m'h' > 



2 
2 
2 



< to,k c > 2 ||/ij| 2 + < /i,k c > 2 ||m|| 2 



:|to'|| 2 ||/i|| 2 + 2 <m,K c xh 2 , mn c > 



|to|| 2 ||ft'|| 2 + 2 < /i,k c >< to 2 , hn c > 



+ ^ < m'*,h 2 > 

+ | <TO 2 ,^' 2 > 



| Hm'/j - to/i'II 2 - ^ (Uto'II 2 + ll^'H 2 ) - |- (< to, « c > 2 + < ft, « c > 2 ) 



4[r 2 (ft,m)-r (to, to) r (ft, ft)] 



— (< to, k c > h+ < h, k c > m) + mhn 
L <P 
i 



— (2 < TO, K c > TO — K c ||TO,|| 2 ^ + TO 2 K C — (2 < h, K c > h - K c \\h\\ 2 ^j + ft 2 K 



(< TO, K c > 2 h 2 + < ft, K c > 2 TO 2 ) 



/ \ 2 



2 < m, k c > rriK c \\h\\ 2 + 2 < ft, k c > hn c \\m\\ 2 — «; 



2 111,112 



TO 
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-21 — 1 < 



to, k c > < ft, k c > mh + — (m 2 K 2 \\h\\ 2 + h 2 K 2 ||to|| 2 ^ 



Taking into account that ft, to form an orthonormal basis, 

Gi (r D (ft, to) , r D (ft, to)) - gi (r (to, to) , r G (ft, ft)) 



l/VV/o 2,_.. .2 IM^ + ^(|| mKc || 2 + 



4 \ 



3 < to, k c > 2 + 3 < ft, k c > z 



4 (f 



||/iK C || 2 ) 



16 



Putting all of this together, we get the sectional curvature 



k c (P(m, h)) = | \\m'h- mh'f - ^ (||m'|| 2 + \\h'f) + ^ (||™k c || 2 + ||/^ c || 2 ) 
V {<m,K c > + <ft,K c > ) - _||/s c || (11) 

Each of the last three terms on the right-hand side is bounded: 

2 1 1 at 1 1 2 2f \\k \\ 2 

\\m Kc \\ 2 + \\h Kc \\ 2 < "^ lo ° , <m,K c > 2 + <h, Kc > 2 < "^ lloc , || Kc || 2 < eW^lt 

where H^dlon = rnaxg |k c (#)|. Therefore, the boundcdness of the sectional 
curvature from above and hence the minimality of a geodesic crucially depend 
on the first two terms. For a fixed m, the magnitude of each of the two terms 
depends on \ \h || which can be made arbitrarily large while keeping \ \h\ \ fixed 
by making h highly wiggly. 



Proposition 7 For a given m, the sectional curvature is bounded from above 
if and only if \\m\\ 2 < or, equivalents, m 2 (9) <.^—m 2 since ^ — ^ ■ 



ip J S ! m 2 d9. 

Proof: We need to estimate only the first two terms on the right-hand 
side of Eq. (11). 

Suppose IjmH^Q < — . Then, 

^\\m'h-mh'f - ^(\\m'f + \\h'f) 



^(m'h) 2 - ^-m' 2 - V mm'hti + ^-( m 2 -^)h' 2 



de 



<f ikiiL-^ / ' y y{h 2 )>de 

<fKH a oo + f ||(mTL<oo 

Conversely, suppose HmH^ > Choose e such that U = {9 : m 2 (9) > 
/f 2 + e} is not empty. Let h be a high frequency wave function with 
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supp(h) C U. Then, 



'rah 



'\\ 2 -^(\\m'f + \\h'f) 



= <pf \^{m'hf - ^fm' 2 - tpmm'hti + ^(m 2 -^) h' 2 
J s i [2 2<p ' 2 \ <P J . 



de 



> 



> 



' 2ip 



||m 



'l|2 



2\//| 



Urn 2 ) 



^\\h>\\ 2 



which tends to oo as the frequency of the wave function h tends to oo. Q.E.D. 
If U = {9 : m 2 (0) > ip' /ip 2 + e} is not empty, 



l = ip 



[ m 2 d6+ [ m 2 d9 

'j^a>(£ + «p)\u\ 



and hence, \U\ < tp/tp'. If ip = £, \U\ < £ and if ip = e M , \U\ < l/A. Thus, the 
case when HmH 2 ^ > ip' /ip 2 may be seen as a generalization of the rectangular 
bump considered in §3. 

If ip = £, the sectional curvature is bounded if and only if m 2 (9) <.m 2 
which is true if and only if m — l/£. Setting m = l/£, we get 



+ 



A£ 2 



<j,K c > 2 + < h,K c > 2 



which is always positive. If h and k c additionally have disjoint supports, the 
sectional curvature equals 3n 2 N 2 /£ 4 where N is the rotation index of C. 

If ip = e M , the sectional curvature is bounded if and only if m 2 {9) <.{A£)m 2 . 
In particular, the sectional curvature is unbounded for every m if £ < l/A. The 
analysis of rectangular bumps in §3 suggests the conjecture that when ip — e Ai , 
a geodesic is locally minimal if and only if a 2 (9) < (A£)a 2 where a = \c^\ (no- 
tation of Eq. 10). 

For an example of a negative sectional curvature, consider the unit square 
with slightly rounded corners. Choose m, h such that supp(m) and supp(h) 
are disjoint and concentrated along the straight portions of the square. Then, 



k c (P(m,h)) = -^-(\\m'\\ 2 + \\h'\\ 2 ) - ^ 



Acknowledgement: Suggestions of David Mumford contributed greatly 
to this paper. 



18 



References 



1. G. Charpiat, O. Faugeras and R. Keriven, "Approximations of shape 
metrics and application to shape warping and empirical shape statistics" , 
To appear in J. of Foundations of Computational Mathematics. 

2. P. Michor and D. Mumford, "Riemannian geometries on spaces of plane 
curves" , |arXiv:math.DG/0312384| v2, Sep. 22, 2004. 

3. E. Klassen, A. Srivastava, W. Mio and S.H. Joshi, "Analysis of planar 
shapes using geodesic paths on shape spaces" , IEEE Trans. PAMI, 26(3), 
pp. 372-383, 2004. 

4. A. Yezzi and A. Mennucci, "Conformal Riemannian metrics in space of 
curves", EUSIPCO04, MIA, 2004. 

5. A. Yezzi and A. Mennucci, "Metrics in the space of curves", 
arXiv:math.DG /0412454| v2, May 25, 2005. 



19 



